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Abstract 

We calculated the SU(2) gluon propagator in Landau gauge on an 
anisotropic coarse lattice with the improved action. The standard 
and the improved scheme are used to fix the gauge in this work. 
Even on the coarse lattice the lattice gluon propagator can be well 
described by a function of the continuous momentum. The effect 
of the improved gauge fixing scheme is found not to be apparent. 
Based on the Marenzoni's model, the mass scale and the anomalous 
dimension are extracted and can be reasonably extrapolated to the 
continuum limit with the values a ~ 0.3 and M ~ 600MeV. We 
also extract the physical anisotropy £ from the gluon propagator 
due to the explicit £ dependence of the gluon propagator. 
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I. Introduction 



The gluon propagator is well defined in the pure non-Abelian gauge theory and is thought 
of one of the most basic quantities of QCD. Due to the asymptotic freedom, the gluon 
propagator at small distances (or large momentum) can be described sufficiently by the 
perturbation theory in which the gluon is expected to propagate as a massless particle. 
However, any asymptotic gluon states and isolated quarks are absent, and this motivated 
the confinement mechanism that gluons and quarks are only confined inside hadrons. This 
physical structure may imply that the gluon propagator is quite different from those of 
stable particles. To understand the structure of the gluon propagator, especially its infrared 
behavior, a nonperturbative study of the gluon propagator is required. 

Lattice simulation provides a powerful method to study the nonperturbative properties of 
QCD. The first lattice simulations of the gluon propagator were carried out by Mandula and 
Ogilvie and by Gupta et al in the late 1980s, in whose studies the gluon propagator in 
Landau gauge was interpreted in terms of a massive particle propagator. The simulations of 
Marenzoni et al gave the result that in terms of an anomalous dimension the propagator 
fell off quite differently from that in perturbation theory. The most recent simulation carried 
out on a very large lattice by Leinweber et al gave the implication that the gluon propagator 
should be interpreted in terms of neither a massless particle nor a free massive one |4]]. 
However, even though there has been substantial progress in the lattice simulation of gluon 
propagator in the last decade, a definitive picture of the gluon propagator has not been 
obtained yet. For reference see the overview of Mandula pfl. 

The development of the improved actions is an important progress in lattice QCD in last 
decade. If the lattice action is properly improved, the costs of the Monte Carlo simulation 
will be greatly reduced and the numerical study can be performed on much coarser lattice 
system. Ma || simulated the SU(3) gluon propagator in Landau gauge on a coarse lattice 
using tadpole improved Symanzik's action with substantially reduced lattice artifact. In his 
work the propagator was fitted using Marenzoni's model and gave the result that the mass 
scale turns out good scaling behavior and the anomalous dimension can be well extrapolated 
to the continuum limit a — ■> 0. 

Since the anisotropic lattice has a lot of advantages in the Monte carlo simulation of 
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lattice QCD on coarse lattice, it is meaningful to study the behavior of gluon propagator 
on anisotropic lattice and the effect of this temporal-spatial asymmetric lattice on the gluon 
propagator because the lattice definition of gluon propagtor is explicitly dependent on the 
anisotropy ratio £. Another goal of this work is to explore the efficiency of the improved guage 
fixing scheme introduced by Bonnet et al |7j . This work is organized as follows. In section II 
we describe the lattice formulations used in this work. Section III is devoted to the simulation 
details where we give the mass scales and the anomalous dimensions extracted from the 
lattice propagator as well as their lattice spacing dependence. In Section IV we discuss the 
effects of the asymmetric lattice on the gloun propagator and extract the anisotropy ratio 
for various lattice spacings. Section V is the summary. 

II. Lattice Formulations 

We use the naive steepest descent method to fix the Landau guage || by maximising the 
usual Landau gauge fixing functional. There are many choices of the definition of the gauge 
fixing functional. The frequently used version (called standard one in this work) is 



f?l{U}} = £ \tt {u^(x) + U°(x?} , (1) 

where 



2' 



U°{x) = G(x)U^x)G(x + /i) f (2) 

and 



G(x) = exp j-i£u; a (>)T a j 



Taking the extremum condition 
6Jf 



Su a (x) 



^r«E uZ(x-a)-u°-(u°(x-ti-u 



A' 



(3) 



A 1 (x) = (4) 



and expanding the gauge links, U^{x) = Vexp \ iag J dtA^(x + fiat) \, in the small lattice 



spacing, one obtains the Landau gauge condition on the lattice as 



(5) 
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The lattice artifact to the gauge condition is of order 0(a 2 ). An alternative gauge fixing 
functional is introduced by Bonnet et al @] as follows, 



Similarly the equation 



5uj a (x) 



« E Tr [ U »( x ~ ^)U^(x -ft)- U^(x)U^x + A) - h.c. 
= A 2 = 



(6) 



(7) 



gives another realization of lattice Landau gauge, 



E<VW = E 



j^x) + 0(a*) 



(8) 



The 0(a 2 ) artifact can be removed by combining the two types of the gauge fixing functional 
to an improved one f7j, T^ mp = \T® — j^-F^i where refers to the tadpole improved 
parameter of the gauge links in \i direction. The improved functional gives the following 
new function 

A 1 

(9) 
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A/ mp (x) = -Ai z - - — A 2 {x). 
3 3Uu 



The steepest descents approach is realized by defining the local gauge transformation as 

a 



Gi(x) 



(10) 



where a is a tuneable step-size prameter. The lattice Landau gauge is fixed by decreasing 
the quantity, 

ei = ^ y ETr{A i {x)A i {*V}, (11) 
to a sufficiently small number. Here the suffix i stands for 1, 2 and Imp. 

In this work, the simulations are performed on an 8 3 x 24 anisotropic lattice by using 
tadpole- improved Symanzik's action ||[10|]. The lattice spacing in the temporal direction is 
smaller than that in the spatial direction with the bare anisotropy ratio £o- The action we 
use is given as 

5 P ss i 1 R xx i 1 iL'., 



= P E 

x,s>s' 

+ /?E 



12 Co< 12 £o«S 
3 u 2 u 2 12 Ugti 2 _ 
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where /3 = g is the gauge coupling. P^ v is the plaquette operator and R^ v 2x1 
rectangular operator, u s and Ut are tadpole improvement parameters. u s is defined by 
u s =< (l/2)TrW sp > 1//4 , while u t is set to be 1 because of u t = 1 — 0(£q 2 ) (£o is much bigger 
than one). 

The continuum gluon propagator D°^ u in Landau gauge is characterized by a single func- 
tion D(q 2 ) as 



D*(q) = 6*\6r-^)Dtf), (13) 

where a and b are the color index. The structure of D(q 2 ) is of physical interest. In the large 
momentum range D(q 2 ) takes its perturbative form and is proportional to 1/q 2 , while its 
infrared behavior may be more complicated and should be reconciled with the confinement. 
The lattice simulation of the gluon propagator is to determine its structure nonperturbatively. 
The lattice version of D(q 2 ) is denoted as Di(q), which can be constructed in terms of the 
Fourier transfromed vector potential A^q). The vector potential of SU(2) gauge field is 
defined from the gauge links 

Ap{x + o M A/2) = ^(U,(x) - Ul(x) - ±Tr(U,(x) - Ul(x))) (14) 
= — A^(x + a M /t/2). 

where is the lattice spacing in fi direction, and its Fourier transformation is defined as 

am = (ik) Y,z~ iq<x+a » m M*+<»/ 2 ) = ( 15 ) 

V /J / x a H 

where q is the momentum 

% = - T n m n M = 0,1,...,L M - 1. (16) 

Here the notation A is adopted for the convenience of the future discussion. The function 
D L (q) is defiend as 

D L {q) = Qr[f 2 r , J2(0\Tr(A,(q)4(l)\°)' ( 17 ) 
which approaches to the function D(q 2 ) in the limit of a M — ► 0. 

III. Lattice Simulation 



We used the action defined in Eqn. ([l2|) to generate gluonic configurations at (3 = 1.0, 1.1, 1.2 
on an 8 3 x 24 anisotropic lattice. The samples of configurations and the values of u s and 

5 



a s [TT|] for various (5 are listed in Table 1. The configurations are separated by 50 Monte 
carlo sweeps, each of which is composed of 4 heatbath and a microcanonical iterations. We 
performed the standard gauge fixing and the improved gauge fixing for each configuration 
respectively and stopped the gauge transfromation at 0i < 1CT 8 . The Landau gauge is surely 
reached for each configuration, for it was observed that for any q = (q, 0), trA^(q)A\(q) are 
of order 10~ 6 and are about five order smaller than other components. 

It is expected that the function D^q) measured on the lattice approximately depend 
only on q 2 = ^lq 2 f° r the small enough a^q^. Previous simulations, which were performed 
on large lattice with the Wilson action, gave the results that this is only true for a small 
range a 2 q 2 < 1. By applying the improved action we find that the function Di(q) from our 
simulation behaves well as a function of q 2 in a much larger momentum range < a 2 q 2 < 9 
which corresponds to almost the full first Brillouin zone (see Figure la, lb, lc, the data 
points are the averages of D^q) with the same q 2 ). This range corresponds to almost the 
full first Brillouin zone. From now on we denote D^q) as D^q 2 ). 

We use the model || to fit the data points of D L (q 2 ), 

Dl{q2) = (irTFF' (18) 

where M is a dimensional parameter, a the anomalous dimension and Z a parameter related 
to g and the renormalization constant of wave-function. We perform the fitting in two ways 
which use the data sets including (fit 1) and excluding (fit 2) the zero momentum point 
respectively, and the results are listed in Table 2 and Table 3. The quality of both fits are 
acceptable due to the very small x 2 /d-o.f. It is found that the effect of the improvement 
of gauge fixing is not apparent, for the fitting results of Di(q 2 ) through the standard and 
the improved gauge fixing procedures are almost the same. The mass scale M in physical 
unit is listed in Table 4. In fit 2, the mass scale is approximately 600MeV within the errors 
and insensitive to the lattice spacing, while in fit 1 the differences of M for various lattice 
spacing are significant. 

The anomalous dimension a is extracted from both the fits. In fit 1 the value of a varies 
with (3 and the a s dependence of a can be well described as 

a(a s ) = a + axaj,, (19) 



6 



with 

a = 0.309(10) ai = -19.9(1.2) X 2 = 0.097. 

In fit 2, the values of a are averaged to be a = 0.282(15) and less related to the lattice 
spacing. 

Even though the small discrepancy, both of the fits give the non-zero anomalous dimen- 
sion a ~ 0.3. This value is compatible with the result of Ma's simulation of SU(3) gluon 
propagator |j)'(a = 0.285(20)) but smaller than the results of Marenzoni || (a ~ 0.4 — 0.6). 

IV. The Renormalized Anisotropy Ratio 

Since the MC simulation is performed on an anisotropic lattice, there naturally arises a 
quesiton: what is the effect of this asymmetry of the temporal and the spatial direction 
on the physical observables and how differently the physical obsevables are affected the 
renormalized anisotropy ratio from the bare parameter £o ( here the bare parameter refers 
to the parameter in the action which justifies that the lattcie action would take the same 
form as that of the continuum action at the classical level when the continuum limit a — ► 
is approached, while the renormalized one, £ = £(£o>/3) is defined by the ratio a s /a t and is 
a funciton of the bare one and the coupling constant owing to the renormalization effects.) 
This question has been explored by Klassen et ai. who extracted £ from the quark static 
potential |I2"|| . We propose here a new scheme to extract the renormalized anisotropy ratio 
from the gluon propagator measured on the anisotropic lattice. 

The lattice version of the funciton Di{q 2 ) and the continuous momentum are explicitly 
related to £ as follows 

D L (q) = * £(Olrr(A M (g)At(g)10) (20) 



oc 



3 



= D L {q,(). 



and 



271 V/V ^ e 



9 2 =h-T + • (21) 



i 
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In practice, the renormalized anisotropy £ is unknown at the beginning of the calculation 
and the measured quantity is Di(q,£,o = 3). For convenience we introduce a notation 
q 2 = H n \ where n M is defined in the Eqn.([l6|). If the function D^(q) is a function of q 2 
and £/£o = 1, the rotational invariance requires that D^(q) take the same value at all the 
momentum with the same q 2 . From the analysis of our results we found that this is only 
correct for the momenuta with the same q 2 and the same n 4 . For those momenta with the 
same q 2 but different n^, the difference of Di(q) is significant. This discrepancy can be 
attributed to the fact that £/£ 7^ 1, as can be found in the following. 

With the conjecture that e/£ is small, only the lowest order term of e = £ — £ is relevant 
in the discussion, D L {q,^Q = 3) is related to the expected D L (q,^) as follows, 

D L (q, = D L (q, 3 - e) = D L (q, 3) - 6e(0|TrA 4 (g)A 4 |0> + 0(e 2 ). (22) 

For simplicity we define 

a(q) = (0\TrA 4 (q)A 4 ( g y\0). (23) 

Another reason for the discrepancy is the £ dependence of the continuous momentum also 
introduces a small shift Aq 2 of momentum square between two momenta q and q' with the 
same q 2 . With both the two factors in consideration and accepting the fact that D^q,^) 
is well described by a function of q 2 , this discrepancy, described by rj = ^^'^ , can be 
re-expressed in the lowest order of e = £, 

D L (q',3) D L (q',Z) + 6ea(q') 



D L (q,3) D L (q,3) + 6ea(q) 



(24) 



DL(q',Q 



1 + 6e 



a(q') a{q) 



D L (q>,0 D L (q,0. 



l-Ag 2 ^lnD L (g 2 ,OI^ 



l + ^J^(S^a(q')-a(q) 
D L (q,0 \D L (q',0 



Since only the lowest order of e is considered, after replacing approximately the ratio D £ ^ ^ 
with I/77 and rearranging the Eqn.(^), the quantity e can be expressed by the quantities 
which can be measured directly from the lattice simulation (see the following equation), 

<q 2 ) = ~, ; — 7 r- (25) 

IS(q 2 )^\nD L (q,0 + ^ (a(q) - \a{q>)) 

" 2 

Here the momenta q and q' satisfy q 2 = (q') with 5{q 2 ) = §e((g') 4 — q\). The results of e(q 2 ) 
for p = 1.0,1.1,1.2 are listed in Table . Theoretically, £/£ should be independent of q 2 . 



From the Table it is found that the measured £/£o are insensitive to q 2 for most momenta. 
However, a few irregular values such as q 2 (S(q 2 ) = 6(1) for (3 = 1.0, q 2 (S(q 2 ) = 4(1) for 
(3 = 1.1 and q 2 {S(q 2 ) = 4(1) for (3 = 1.2, are a bit larger than the others. We attribute 
it to the reason that the rotational invariance is violated to some extent so that the value 
of Di(q) with momenta directed along one of the four axes is different from those with 
momenta directed off-axis. Averaging the regular values we obtain 

£/£o = 1.23 for (3 = 1.0 

£/£o = 1.15 for (3 = 1.1 

£/£ Q = 1.16 for (3=1.2. 

These are raw results at the lowest order of e and we do not plan to extract the (3 dependence 
of £/£o from only three data points. What we can conclude at present is that the renormal- 
ization effects for the anisotropy ratio is significant and our results is consistent with those 



from the study of static quark potential by Klassen et al |12]. Morningstar et al [13[] have 
simulated the SU(3) glueball mass on the anisotropic lattice with the same action as that we 
use and reported that £/£o extracted from the static quark potential is apporximately one 
within few percents discrepancy. Work is underway to explain the reason why the anisotropy 
ratios extraceted from different quantities do not comply with each other. 

V. Summary 

We studied the SU(2) gluon propagator in Landau gauge on an 8 3 x 24 anisotropic lattice 
by using the tadpole improved Symanzik's action. Although the improvement scheme is also 
used to fix the gauge, we found that its effect is not significant, while the lattice artifact is 
significantly suppressed by using the improved action. The function D^q) in the propagator 
can be well described as a function of continuous momentum square q 2 even in the full first 
Brillouin zone. We used the model suggested by Marenzoni et al to fit the function Di(q) and 
extracted the mass scale and the anomalous dimension therein. It is found that the value 
of the mass scale is approximately independent of the lattice spacing and the anomalous 
dimension can be well extrapolated to the value about 0.3 in the continuum limit. 

The effects of the asymmetric lattice on the gluon propagator are also explored in this 
work. Based on the fact that the lattice definition of the gluon propagator is related to the 
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renormalized anisotropy ratio £ = a s /a t , we deduced a formula to measure the £ nonpertur- 
batively at various momentum. Even in the lowest order of e = £o — £> the measured values 
of £ are insensitive to the momentum and are consistent with the results extracted from the 
static quark potential. 
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Tables and Figures 
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p 


U s 


a s (fm) 


No. of Configs. 


1.0 


0.8215 


0.277 


200 


1.1 


0.8365 


0.233 


200 


1.2 


0.8593 


0.183 


200 



Table 1: The parameters used in the simulations of this work. 





Standard fixing 


Im 


Droved fixing 


p 


ajM 2 


a 


X 2 /dof 


a 2 s M 2 


a 


X 2 /dof 


1.0 


0.552(24) 


0.191(15) 


0.17 


0.547(24) 


0.203(15) 


0.17 


1.1 


0.357(14) 


0.254(13) 


0.17 


0.360(14) 


0.268(13) 


0.15 


1.2 


0.311(12) 


0.285(13) 


0.13 


0.313(11) 


0.298(13) 


0.12 



Table 2: The results of the fit 1 which includes the zero-momentum data point, a and a 2 s M 2 
in the propagator model are extracted at various value of f3. The data illustrate that the 
effect of the improved gauge fixing scheme is not apparent. 





Standard fixing 


Im 


Droved fixing 


p 


a 2 M 2 


a 


X 2 /dof 


a 2 s M 2 


a 


X 2 /dof 


1.0 


0.808(95) 


0.273(33) 


0.11 


0.807(88) 


0.288(32) 


0.11 


1.1 


0.454(70) 


0.289(28) 


0.18 


0.439(69) 


0.297(28) 


0.17 


1.2 


0.309(66) 


0.284(27) 


0.15 


0.317(66) 


0.299(27) 


0.13 



Table 3: The results of the fit 1 which excludes the zero- momentum data point, a and a 2 M 2 
in the propagator model are extracted at various value of /3. The data illustrate that the 
effect of the improved gauge fixing scheme is not apparent. 
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p 


1.0 


1.1 


1.2 


M(MeV)(Fit 1) 


529(12) 


506(10) 


601(12) 


M(MeV)(Fit 2) 


640(37) 


571(44) 


599(64) 



Table 4: Tha mass scales in physical unit. In fit 2, within the range of errors, the mass scale 
takes approximately the same value and is independent of the lattice spacing. 





p = 1.0 


P = 1.1 


P = 1.2 




e 


£/£o 


e 


£/£o 


e 


£/£o 


1(1) 


-1.04 


1.35 


-0.31 


1.10 


-0.40 


1.13 


2(1) 


-0.69 


1.23 


-0.52 


1.17 


-0.50 


1.17 


3(1) 


-0.79 


1.26 


-0.32 


1.11 


-0.27 


1.09 


4(1) 


-0.66 


1.22 


-1.14 


1.38 


-1.09 


1.36 


4(4) 


-0.54 


1.18 


-0.40 


1.13 


-0.39 


1.13 


5(1) 


-0.75 


1.25 


-0.40 


1.13 


-0.50 


1.17 


5(4) 


-0.63 


1.21 


-0.41 


1.14 


-0.43 


1.14 


6(1) 


-1.09 


1.36 


-0.53 


1.18 


-0.58 


1.19 


6(4) 


-0.68 


1.23 


-0.50 


1.17 


-0.49 


1.16 


7(1) 


-0.85 


1.28 


-0.44 


1.15 


-0.62 


1.21 


8(4) 


-0.71 


1.24 


-0.63 


1.21 


-0.60 


1.20 



Table 5: The difference of the renormalized anisotropy ratio and the bare one e = £o — £ 
is measured at various q 2 . The ratio £/£o is insensitive to the momentum except for a few 
irregular values. These irregular values may be the result of the significant violation of the 
rotational invariant. 
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Figure 1: The lattice gluon propagator D^q) versus the dimensionless continuous momen- 
tum a 2 q 2 measured at various (3. The filled circles are the data points from the simulation, 
the dash line shows the result of fit 1 and the full line the results of fit 2. a) (3 — 1.0; b) 
(3 = l.l; c) 13 = 1.2. 
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0.02 0.04 0.06 0.08 0.1 

a 2 



Figure 2: The lattice spacing dependence of the anomalous dimension a is illustrate in the 
plot, a through fit 1 (filled circles) can be extrapolated to the continuum limit by the 
model a(a) = 0.309 - 19.9a 4 with the x 2 = 0.097 while a through fit 2 (filled square) is 
approximately independent of the lattice spacing within the errors. 
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